Intrinsic susceptibility and bond defects in the novel 2D frustrated antiferromagnet 

Ba2Sn2ZnCr7pGaio_7p022 
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We present microscopic and macroscopic magnetic properties of the highly frustrated antiferro- 
magnet Ba2Sn2ZnCr7pGaio-7p022, respectively probed with NMR and SQUID experiments. The 
T-variation of the intrinsic susceptibility of the Cr'^"'" frustrated kagome bilayer, Xkag, displays a max- 
imum around 45 K. The dilution of the magnetic lattice has been studied in detail for 0.29 < p < 0.97. 
Novel dilution independent defects, likely related with magnetic bond disorder, are evidenced and 
discussed. We compare our results to SrCr9pGai2-9pOi9. Both bond defects and spin vacancies do 
not affect the average susceptibility of the kagome bilayers. 

PACS numbers: 75.30.Cr, 75.50.Lk, 76.60.-k 
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Geometric frustration in magnetism has proven to 
yield various original ground states in the past decade, 
from RVB "spin hquids" , to exotic freezings as spin ice. 
Most of these systems have in common a simple Heisen- 
berg Hamiltonian with near neighbor (nn) antiferromag- 
netic interactions on a corner sharing lattice, e.g. the 
kagome lattice in two dimensions Q. In the latter case 
and S = ^, the ground state is expected to be the long 
sought resonating valence bond (RVB) state 0- This 
opens not only the way to novel states in low dimensional 
quantum magnetism [3| but the RVB state has been also 
advocated in manyfields of condensed matter such as 
high Tc cuprates |j| or recently discovered cobaltites, 
where frustration is inherent to the triangular Co net- 
work j^]. Geometric frustration is also a rapidly devel- 
oping field since in the classical Heisenberg case, small 
perturbations such as disorder, anisotropy or weak ad- 
ditional interactions, lift the frustration-induced macro- 
scopic degeneracy of the ground state . 

So far, there is no ideal experimental candidate for a 
spin liquid behavior, but a few, like the newly discovered 
Ba2Sn2ZnCr7pGaio-7p022 compound (BSZCGO(p)) 0, 
retain the essence of the expected features, e.g. a short 
correlation length Q which prevents the occurence of 
any long range order in a broad T-range below the 
Curie- Weiss (CW) temperature 9cw, a large residual 
entropy j3|, and persistent fluctuations at T 0- 
The corner sharing magnetic lattice of BSZCGO(p) is 
made of well decoupled Cr^+ {S = |) 2D kagome bilay- 
ers only. Interestingly, the nn direct couplings between 
Cr'^+ ions make the Hamiltonian quite close to the ideal 
Heisenberg case ^ spin-glass like transition occurs 
at Tg « 1.5 K 0, well below Ocw ~ 350 K, which yields 
the highest frustration ratio f = dew /Tg ^ 230 reported 
so far in a compound where frustration is driven by cor- 
ner sharing equilateral triangles Q. 

In comparison, the long-studied parent Cr'^+ com- 
pound SrCr9pGai2-9pOi9 (SCGO(p)) [ill (/ ^ 150) is 
not made of frustrated bonds only, since | of the Cr^+ 



form pairs which separate the bilayers and couple into 
singlets below 216 K Ga^+ /Ct^+ substitutions of 

the order of a few percent {= I — p) represent the major 
drawback in SCGO(p) since they introduce paramagnetic 
centers in between the bilayers by breaking these Cr'^+ 
pairs. This not well-controlled parameter, together with 
the shorter inter-bilayer distance, 6.4 A instead of 
9.4 A in BSZCGO(p) 0, is hkely to make the physics 
of SCGO(p) less 2D and complicates the search for the 
intrinsic kagome physics. 

We present the first NMR study of BSZCGO(p) and 
show that the intrinsic susceptibility (xfcag) displays a 
maximum at a temperature comparable to the nn Heisen- 
berg coupling J. The similarity of the T-variation of Xkag 
with the case of SCGO(p) stresses the universal character 
of this maximum in Xkag for kagome bilayers. As a local 
probe, Ga NMR, contrary to SQUID, offers a unique way 
to disentangle Xkag from the contribution of "magnetic 
defects" to macroscopic susceptibility 0, 0| . This gives 
evidence for two types of defects, for 0.29 < p < 0.97, 
(i) the spin vacancies due to the Ga'^"'"/Cr^+ substitu- 
tion which was shown to generate a staggered response of 
the surrounding magnetic background [l0| and (ii) novel 
dilution independent defects, discussed here. 

All the samples were synthesized by a solid state reac- 
tion of BaCOs, Cr203, Ga203, ZnO and Sn02 in air at 
1400°C. They were characterized by X-ray diffraction, 
macroscopic susceptibility (xmacro) measurements and 
^^Ga NMR. The synthesis of pure samples with p > 0.97 
failed. Indeed, in this case. X-ray diffraction shows a 
Cr203 parasitical phase and an extra kink of Xmacro is 
also observed around 6 K. This sets the maximum of the 
Cr'^+ lattice occupancy to a value of Pmax = 0.97. Below 
Pmaxi Ga substitutes randomly on the Cr(la, 6i) site. 

The two Ga'^+ site environments {p = 1), are displayed 
in the inset of Fig. ^b). The Ga(2d) site is coupled to 3p 
Cr(la) and 9p Cr(6i) through neighboring oxygens. Its 
local charge environment is approximately cubic, whereas 
the Ga(2c) site belongs to a tetrahedron elongated along 
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FIG. 1: ^^Ga spectra, a) T > 80 K. The ''^Ga(2c) first order 
quadrupolar contribution appears as a flat background in this 
fleld range, b) T < 40 K, lines are broadened, o and □ are 
the short r and the rescaled long r spectra respectively. A is 
for the Ga(2d) contribution, given by their subtraction. The 
dotted arrows point the center of the Ga(2d) line. Continuous 
(dashed) lines are gaussian broadened ''^Ga{2d) C^Ga{2c)) 
powder pattern quadrupolar simulation, with uq = 3.5 MHz 
(12 MHz) and rj = 0.6 (?? = 0.04). Inset : Ga^+ sites. 



the c-axis and couples to 3p Cr(6i). 

We performed ^^Ga sweep field NMR experiments at 
constant frequency i^i = 84.365 MHz, using a f -t-tt spin 
echo sequence. Two sites can be clearly identified in 
Fig. n with the expected ratio of integrated intensities 
close to 2:1. At 150 K, the double peak around the 
reference field Hi « 6.498 T in Fig ^a) is character- 
istic of a 0.53% slightly shifted 5^—5 nuclear tran- 
sition, with second order quadrupolar powder pattern. 
The corresponding first order satellite singularities are 
clearly seen in the less expanded spectra of Fig. ^b). 
From these singularities, we extract a quadrupolar fre- 
quency ^^t'Q ~ 12 MHz and an asymmetry parameter 
77 = 0.04(2) Q. The other hue is more shifted ~ 2% 
(Fig. ^a)). Its shape can be explained by a broadened 
5 = 1 quadrupolar pattern, with ^^i^q = 3.5 ± 1 MHz 
(Fig. mb)) [251. The intensity ratio, the very different 
quadrupolar frequency associated with different charge 
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FIG. 2: Shift of the ''^Ga{2d) (Ga(4f)) line for BSZCGO 
(SCGO). Inset : full scale down to 0. Lines are fits (see text). 

environment and the analogy with SCGO allow us to 
identify unambiguously the Ga(2c) and Ga(2d) sites. 
Though, contrary to SCGO(p), we could not identify any 
line typical of substituted Ga. 

At low-T, the two lines are found to overlap due to 
strong broadenings. In order to deconvolute the two con- 
tributions, we took advantage of the different transverse 
relaxation times T2 of the two lines (Fig. ^). The first 
step was to use a large r spectrum (r = 200 ^s) to isolate 
the Ga(2c) contribution (T'j2c ~ 200 ^s). Its linewidth 
AH is extracted using a gaussian convolution of the high- 
T quadrupolar pattern. In the second step, we used the 
short T spectrum (r ^ 10 fis) where the Ga(2d) contribu- 
tion is evident, and subtracted the Tjjd normalized long 
T spectrum. The width and the shift of the Ga(2d) line 
can then be extracted from a fit using a ''^j^q = 3.5 MHz 
quadrupolar pattern convoluted by a gaussian line shape, 
not affected by the 30% uncertainty on "^^vq. As a mat- 
ter of consistency, the ratio /J^ /^J^ = 2 ± 0.8 stays close 
to 2/1 down to 10 K. No significant loss of intensity was 
observed down to this temperature. 

The kagome bilayer susceptibility Xkag is revealed in 
the shift Kl^ of the Ga(2d) line [ij , presented in Fig. 
It increases when T decreases for T > 80 K, reaches a 
maximum around T^ax ~ 45 K and decreases below. A 
phenomenological CW law, K^^ — Cj^mrI^ + &nmr), 
is an accurate fit for T > 80 K with Onmr = 380 ± 10 K. 
Although very common, this is a rather crude approx- 
imation in a T-range of the order of Ocw- In the ab- 
sence of any prediction for S = |, we used the high-T 
series expansion of Xkag, derived for a S = ^ kagome 
lattice, to correct 9nmr by a factor 1.5 [l^- From 
(^cw = zS{S + 1) J/3 and z = 5.14 the average number 
of nn for the Cr^+ of the bilayer, we extract J « 40 K. 

In SCGO(0.95) [l^, a value 0nmr = 440 ± 5 K [ll, 
was found and Cnmr is 20% larger. The small variation 
in the CW behaviors can be explained by a slightly larger 
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FIG. 3: Magnetic contribution to the NMR linewidth 
AH/HiC^Ga{2d)) (•). Data for the ^^Ga(2c) (□) arerescaled 
by a factor 6, corresponding to the ratio of the couphng con- 
stants, consistent with the high-T shifts. The lines are oc 1/T 
fits. Inset : p-dependence of A_ffJJ (AHlj) at 50 K for 
BSZCGO(p) (SCGO(p)). The lines are guides to the eye. 



C]-3+_(^j.3+ coupling, due to a 0.01 A mean shorter bond, 
and A J/ Ad ~ 450 K/A [T^. 

The existence of a gap A in Xkag was already an impor- 
tant issue for SCGO(0.89) [3, since it would be a unique 
signature of the existence of a singlet ground state. For 
S" = ^, its value is predicted to be A ~ J/20 Il6l. which 
should be little influenced by the spin value [23. Its 
definite observation would then require T < A in or- 
der to observe the exponential decrease of Xkag- Since 
our data is masured above J/4 ~ 10 K, we only ten- 
tatively fitted with the phenomenological function 
cx exp(-A/T)/(r + 9nmr)- Although the value A - 
J/10 ~ 4 K is emerging, further details of the variation 
of the susceptibility are not accounted for, e.g. the rela- 
tive sharpness of the variation of Xkag around T^ax ■ 

Rather than a sign of a gap, it was suggested in 
[1^ that the maximum of Xkag would be the land- 
mark of a moderate increase of magnetic correlations, 
later confirmed by neutron measurements [l^ . Suscep- 
tibility calculations were further performed in kagome 
and pyrochlore lattices with Heisenberg spins, using 
the so-called "generalized constant coupling" method 
(GCC) [13. The fits for the 5 = § kagome and py- 
rochlore lattices, in between which the kagome bilayer 
susceptibility is expected to lie, yield very close values of 
the CW temperature Ocw = 238 ± 3 K (Fig. EJ. One 
extracts J « 37 K for BSZCGO(0.97), consistent with 
the high-T series expansion. 

As expected, the GCC computation fits quite well the 
data around T^ax- It only shows that the strongest un- 
derlying assumption of this cluster mean-field approach 
is relevant, i.e. the spin-spin correlation length is of the 
order of the lattice parameter. This is also proven by 
neutron diffraction [2(j|. 

The linewidth is the second important piece of infor- 



mation which we extract from our data. It is indeed 
due to a distribution of susceptibility which is related 
to the existence of defects and to their nature In 
Fig. 13 we compare the linewidth AH /Hi obtained for 
BSZCGO(0.97) to the one obtained for a comparable di- 
lution in SCGO(0.95). In both samples, AH/ Hi follows 
a Curie-like behavior which was shown to be character- 
istic of magnetic defects in SCGO(p). Surprisingly, in 
BSZCGO(0.97), the linewidth is ^ times larger than in 
SCGO(0.95). This ratio cannot be explained by a vari- 
ation of hyperfine constant, which, from the ratios of 
Cnmr, was estimated to a maximum of 20%. In addi- 
tion, the Ga(2d) linewidth in BSZCGO(p) could be mea- 
sured for various p at T = 50 K (Fig.Ofinset)). Contrary 
to SCGO(p), where AH /Hi extrapolates to for a per- 
fect lattice {p = 1), wc find that it reaches an asymptotic 
non-zero value in BSZCGO(p) for 0.86 < p < 1. One 
could wonder whether the dilution of the magnetic lat- 
tice is larger than the nominal one. This can be ruled 
out since (?) we found the expected evolution of the line- 
shape at 300 K with p jij] {ii) muon spin relaxation mea- 
surements indicate a regular evolution of the dynamical 
properties with p 0. We then conclude that the extra 
width is not related to spin vacancies and reflects new 
p-independent defects in BSZCGO(p). 

Defect terms could also be probed using low-T macro- 
scopic susceptibility measurements. We performed mag- 
netization measurements under a field of 100 G, for 
0.29 < p < 0.97, down to 1.8 K. No difference was 
observed between Field Cooled and Zero Field Cooled 
branches down to this temperature, in agreement with 
Along the common route, one can simply analyze 
the data using the phenomenological fit : 



Xr. 



Ccw{p) , Cc{p) 



T + e{p) T + 0c{p) 



(1) 



The first term reflects the high-T weak variation of the 
susceptibility due to AF interactions. From this fit, we 
can extract a p-independent effective moment per Cr^+ , 
Peff — 4.1 ± 0.2//XB, close to 3.87/fiB expected for 
5 = |. From the linear variation of 9{p) ~ 1.5 9cw{p) 
(Fig. ^a)), we extract J ~ 40 K, consistent with our 
NMR result. 

The second term in Eq.^is introduced to fit the low-T 
behavior associated with the defects 0, . It is close 
to pure Curie law, with 9c{p) = 1 ± 0.5 K. Values of 
the Curie constants Cc are reported in Fig. ^b). We 
find a p-independence qualitatively similar to the NMR 
width data and a finite value for p — > 1. To give an 
order of magnitude, this term represents 20% of 5 = ^ 
paramagnetic spins with respect to the number of Cr''+ 
spins. Once the dilution defect term becomes significant 
(p < 0.8), the Cn{p) slope is qualitatively comparable in 
both systems j28j . 

To summarize, low-T Xmacro and NMR width data 
can be both satisfactorily accounted for, only if a novel 



4 




P P 



we showed that bond disorder is likely to give raise, like 
spin vacancies, to extended defects which underline the 
correlated nature of the kagome network. This study of 
"impurities" might be a major avenue for revealing the 
nature of the ground state in a similar manner as in low- 
D quantum magnetism studies. 

We thank F. Bert, J. Bobroff, A. Keren, C. Lhuillier, 
R. Moessner, C. Payen and P. Sindzingre for discussions. 



FIG. 4: Fits of the SQUID data (Eq.p. a) Linear fit of 9{p) 
for BSZCGO(p). b) Cc{p), lines are guide to the eye. 

p-independent defect-like contribution is taken into con- 
sideration, a result somehow surprising in view of the 
close similarity between kagome bilayers. We now elab- 
orate about the nature of these defects. The only major 
change lies in the 1:1 random occupancy of the 2d site by 
Ga-^^ or Zn^+ ions which induces different electrostatic 
interactions with the neighboring ions. As an example, 
distances to in a tetrahedral environment vary from 
i~Ga^+-o'^- = 0.47 A to r^„2+_Q2- = 0.60 A which is at 
the origin of a change of the average Ga(4f,2d)-0 bonds, 
from 1.871 A for SCGO |l3 to 1.925 A in BSZCGO 0. 
Similarly, one expects that the Cr'^+ will be less repelled 
by Zn^+ than by Ga'^+ , which certainly induces magnetic 
bond-disorder, i.e. a modulation of exchange interac- 
tions J between neighboring Cr'^^. As seen before, one 
can estimate that a modulation as low as 0.01 A in the 
Cr^+-Cr'^+ distances would yield a 10% modulation of J. 

The presence of non-perfect equilateral triangles clas- 
sically induces a paramagnetic cornponent in the sus- 
ceptibility of the frustrated units j23|. The fact that 
we do not observe any extra Curie variation as com- 
pared to SCGO(p) in the average susceptibility, probed 
through the NMR shift K^^ (Fig.Hl), indicates that such 
unbalanced exchange interactions only induce a stag- 
gered response in the same manner as spin vacancies . 
Whether this could be connected with localization of sin- 
glets in the vicinity of defects and a staggered cloud 
around should be more deeply explored. This might high- 
light the relevance of a RVB approach to the physics of 
the kagome network. Further insight into the exact na- 
ture of the defects likely requires a better determination 
of the bond disorder through structural studies at low-T, 
to avoid the usual thermal fluctuations at room T. 

To conclude, one of our major findings is the maximum 
of the local susceptibility Xkag observed around 45 K for 
both BSZCGO(p) and SCGO(p), which appears very ro- 
bust to the presence of defects either spin vacancies as 
observed in \m or bond disorder as pointed out above. 
This intrinsic susceptibility is very close to a model of 
S' = I kagome or pyrochlore with correlation length of 
the order of magnitude of the Cr'^"''-Cr'^+ bond. Finally, 
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(Fig. 121 inset)). We attribute this difFerence to the extra sured in Xmacro and not in Ti'JJ oc Xkag (unpublished). 

Ga-Cr substitutions in the Cr(4fi,i) sites in SCGO, mea- 



